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This is a set of introductory lectures on the behaviour of a directed polymer in a random medium. 
Both the intuitive picture that helps in developing an understanding and systematic approaches for 
quantitative studies are discussed. 
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I. DIRECTED POLYMERS 

A physicist's polymer is a long string like object, de- 
void of any specific inner structure. It can be treated as a 
random walk on a lattice or as an elastic string in contin- 
uum. The lattice problem is well-defined though in many 
situations the continuum version is easier to treat. For 
the lattice problem, the sites visited are the monomers 
with the steps as bonds. The monomers are indexed by a 
coordinate z which becomes a continuous variable in the 
continuum model. See Fig. ^ 

An important quantity for a polymer is its size or the 
spatial extent as the length N becomes large. For a trans- 
lationally invariant system with one end (z = 0) fixed at 
origin, the average position at z = iV is zero but the size 
is given by the rms value 



(1) 



with ly = 1/2, for the free case. In general the size ex- 
ponent defines the polymer universality class and it 
depends only on a few basic elements of the polymer. In 
addition to this geometric property, the usual thermody- 
namic quantities, e.g. free energy (or energy at T = 0), 
are also important, especially if one wants to study phase 
transitions. 




FIG. 1: (a) A random walk in d dimensions with z as the vari- 
able along the contour of the polymer i.e. giving the location 
of the monomers, (b) Directed polymer on a square lattice. 
A polymer as of (a) can be drawn in d-|- 1 dimensions. This is 
like a path of a quantum particle in nonrelativistic quantum 
mechanics. 

We consider the problem of a polymer where each 
monomer sees a different, independent, identically dis- 
tributed random potentials. Geometrically this can be 
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achieved if the monomers hve in separate spaces. One 
way to get that is to consider the polymer to be a d + 1 
dimensional string with the monomers in d dimensional 
planes but connected together in the extra dimension. As 
shown in Fig. ^ this is a polymer which is directed in one 
particular direction. Hence the name directed polymer. 

For a directed polymer, the size would now refer to the 
size in the transverse d-directions and so Eq. refers to 
the transverse size as the length in the special z-direction 
increases 

A. Why bother? 

The major issue for the random problem is to see if 
the universality class of the pure polymer changes to a 
new one in presence of disorder and if any new phase 
transition can take place. 

Because of randomness (or disorder), any quantity of 
interest has to be a stochastic variable (realization de- 
pendent). Therefore, an additional disorder averaging 
needs to be done over and above the usual thermal av- 
eraging for each realization. For thermal fluctuations, 
there is generally no need to go beyond second moments 
or cumulants (related to various response functions like 
the specific heat, susceptibility etc.) but that cannot be 
said offhand for the disorder averaging. The random- 
ness is put in by hand and is not thermalizable. There 
are possibilities of rare events for which it is necessary 
to distinguish between the average value and the typical 
(i.e. most probable) value. In such situations higher mo- 
ments of the quantity concerned become important. It is 
this aspect that makes disorder problems interesting and 
difficult. 

^^^^^^^^ Example: Typical vs. average 
Consider a random variable x that takes two values 

Xi = e"^ and X2 = e''^,/? > 1, (2a) 

with probabilities 
pi = 1 - e"^, and p2 = e"^. (2b) 

In the limit N 00, the average value [x]^^ — > e^''^-'^^^ 
while the typical or most probable value is a; = Xi with 
probability 1. On the other hand [lna:;]JJ.^, — s- a\/N in 
the same limit, showing that [Inxjjj^ is determined by 
the typical value of the variable while the moments are 
controlled by the rare events. Note that this peculiarity 
disappears if x has a smooth probability distribution in 
the sense of no special or rare events. 

The significance of directed polymer lies in the fact 
that the pure system is very well understood and exactly 
solvable in all dimensions. Polymers show critical- like 
behaviour (e.g., power laws for size) without any need 
of any fine tuning. One hopes that the understanding 
gained in this topic on the role of disorder may be useful 
in other situations as well. Wistful sighs... 

Background: Familiarity with the basic ideas of 
renormalization group or critical phenomena (Ref. 0) 
and polymers (Refs. QjH) would be useful. 



II. HAMILTONIAN AND RANDOMNESS 

A. Pure case 

Taking the polymer as an elastic string, one may define 
a Hamiltonian 

which gives a normalized probability distribution of the 
position vector r at length z 

^^''^)=(2^^"'^''" (i^rf/fcBr = i). (4) 

For the lattice random walk, there is no "energy" and the 
elastic Hamiltonian of Eq. |31just simulates the entropic 
effect at non-zero temperatures. One needs to look at 
the lattice problem in case one is interested in low or 
zero temperature bahaviour. 

For a polymer of length the probability distribution 
gives 

(r)=0, {r^)=N, {Kd/kBT = l) (5) 

so that the transverse size of the polymer is given by 

i?o - (r2)i/2 - AT'' with 1/ = 1/2. (6) 

The power law growth of the size of a polymer as the 
length increases in a reflection of the absence of any 
"length scale" in the Hamiltonian. 

Exercise II. 1 Put in d, K and T in Eqs. 
Get Eq. \^from Eq. 

B. Let's put randomness 

1. Hamiltonian 

Let US now put this polymer in a random medium. In 
the lattice model of Fig. each site has an independent 
random energy and the total energy of the lattice poly- 
mer is the sum of the energies of the sites visited. In 
continuum, the Hamiltonian can be written as 

H = Ho+ dz r]{r, z) S{r{z) - r) (7) 
Jo 

where 77(r, z) is an identical, independent gaussian dis- 
tributed random variable with zero mean and variance 
A > 0, 

b7(r, z)],, = 0, [77(r, z) r^{v\ z')]av = A5(r - v')5{z - z'). 

. (8) 

With this distribution of random energies, we see 
= Hq and so the average Hamiltonian is not of 
much use. 
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2. Partition function 

The partition function for a polymer in a random 
medium or potential is given by 

Z = j Vn e^'^". (9) 

This is a symbolic notation to denote sum over all con- 
figurations and is better treated as a continuum limit of 
a well-defined lattice partition function 

Z^Y. e"''"*'"'^^ (10) 

paths 

where the sum is over all possible paths of N steps start- 
ing from r = at z = 0. and Zq is the partition function 
of the free walker. It often helps to define the partition 
function such that Z{{r] = 0}) = 1 to avoid problems of 
going to the continuum limit (see Eq. ^J. This is done 
by dividing (or normalizing) Z by Zo = /-t^, being 
the partition function of the free walker with /i as the 
connectivity constant (= 2 for Fig. ^p). 

3. Averages 

Because of the randomness, there is no unique par- 
tition function but rather a probability distribution of 
the partition function, P{Z), is needed. Any quantity 
of interest needs to be averaged over such a distribution. 
This averaging is to be called sample averaging, denoted 
by [■•■lav (^s opposed to thermal averaging, denoted by 
(...)). Consequently, sample to sample fluctuations may 
become important and would be a new class of quanti- 
ties which are not meaningful in the pure thermal case. 
The extreme situation is the possibility of "rare" events 
which require distinctions between typical and average 
behaviour. 



C. What to look for? 

For the pure case (A = 0), there is no "energy", only 
configurational entropy of the polymer. But with A ^ 0, 
there may be one or more lowest energy states. What is 
the nature of the ground state? 

Exercise II. 2 For the lattice problem, the energy is the 
sum of iV(— > oo) random energies of the sites visited. 
This is tantalizingly similar to what one needs for apply- 
ing the central limit theorem. Is it so? 

D. Size: disorder and thermal correlation 

So far as the geometrical properties are concerned, we 
first note the lack of translational invariance for a partic- 
ular realization of disorder and therefore {bfr) ^ 0, but 



on averaging over randomness, translational invariance 
will be restored statistically and so [{bfr)],^.^ has to be 
zero. One may therefore consider the size of the polymer 
by the "disorder" correlations 

Cdis = [(&/r')]av or Cdis = [|(6/r)|]L with Qis ~ N^" 

(11) 

Furthermore, the usual correlation function is the ther- 
mal correlation 

Ct ^ [(fo/r')],. - mr)\,. (12) 

These are defined at all temperatures. Of course, the 
probability distribution P{{bfr)) is also of importance. 



1. Free energy fluctuation 

Apart from these geometric properties one would like 
to know the fluctuations in free energy at any non-zero 
temperature (or ground-state energy at T = 0). The 
probability distributions for these quantities would be 
the ideal things to look for. In absence of knowledge of 
this distribution, one may look for various moments and 
fluctuations. 

Let us define F{r, N) ~ — T In Z{r, N), a restricted free 
energy that the end at z = is at r. This free energy, 
being realization dependent, is also a stochastic quantity. 
We may then define a correlation function 

Cpir, N) = [(F(r + ro, N) - F(ro, N))\,, (13) 

which may be assumed to satisfy a scaling form 

CFir,N) = N'U[^). (14) 

Note that the exponent 6 has no meaning in the pure case 
where 6 — 0. This is a new quantity for the disordered 
system. P{r,z) of Eq. (@J, if used for Z{r,z), satisfies 
Eq. (O with 6* = and 1/ = 1/2. 

For r <C A^"^, Cp should be a function of r. Taking 
f(x) ~ x'^ for X ^ 0, we need p ~ Q jv and get 

CF{r:,N) = r^"/", for r < TV' 

^ N^^. for r > N". (15) 

2. Our Aim 

Our aim is two-fold: to get the quantitative estimates 
of v and 9, and to develop an intuitive picture. In fact, 
there are quantitative changes in both the exponents but, 
more importantly, those changes are consequences of a 
drastic change in the qualitative behaviour of the poly- 
mers. At the end, however, a simple picture borrowed 
from pure polymers remain valid. 

On Notation: 
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• To avoid proliferation of symbols, we reserve the 
symbol f to denote an arbitrary or unspecified func- 
tion, not necessarily same everywhere. 

• Sample averaging is denoted by [■■■]av while thermal 
averaging is denoted by (...). 

• The Boltzmann constant is set, most often, to one, 
kB = 1. 

• "Disorder" and "randomness" will be used inter- 
changeably. 



3. Digression: Self-averaging 

Our interest is in large systems. A question of building 
up such a large system arises. In the pure case, we may add 
blocks systematically (see Fig. to generate a unique way of 
going to the limit. For a random system, it is not so straight- 
forward. For a given block A, there are many possibilities for 
each of B,C,D..., well, calling for a probabilistic description. 
Even the issue of thermodynamic limit has to be probabilistic 
in nature; we may expect In Z to be extensive for large size 
with probability one but not necessarily for each and every 
realization. A quantity of import would be the probability 
distribution (sample to sample variation) of the free energy 
itself, P{\nZ). If P(lnZ) is sharply peaked for large size, 
then any large sample would show the average behaviour. A 
quantity with this property is often called self-averaging. This 
may not be the case if the distribution is broad especially in 
the sense discussed in Sec. 1.1. A self-averaging quantity has 
the advantage that one may study one realization of a large 
enough system without any need of further disorder or sample 
averaging. 












FIG. 2: Adding blocks to build a larger system. 



This is equivalent to the statement that the sum over a large 
number of subsystems gives the average value, something that 
would be expected in case the central limit theorem (CLT) is 
applicable. But for many critical systems CLT may not be 
applicable. 

A more practical procedure for testing self- averaging be- 
haviour of a quantity X is to study the fluctuations a|f = 
[A^^lav ~ [^lav smd then check if 

2 

Rx,N = 7^;^ -> 0, as iV ^ oo. (18) 

One may define more general measures like Rx^'n ~ 
a"^V[^]av where cr^' is the pth cumulant. A quantity is 
not self-averaging if the corresponding R does not decay to 
zero. A very slow decay ("weakly" self averaging) would also 
signal practical difficulties. A thumb rule would be that an 
extensive quantity is generally self- averaging, others may not 
be. 



III. Q.: IS DISORDER RELEVANT? 

Let us first sec if disorder is at all relevant. 



A. Annealed average: low temperature problem 

The averaged Hamiltonian is of no use but what about 
the average partition function [Z]^^7 This is called an 
annealed average. With a gaussian distribution for the 
random energy, from Eq. H10() . 

[Zl, = exp(/32AiV/2)cxp(iVln^), (19) 

so that F/N = -T[\nii -\- The entropy obtained 

from this partition function (5 — ~dF / dT) by definition 

has to be positive. This gives a limit T >Tp^ = \J 2^^ - 
Annealed averaging will not work at very low tempera- 
tures. 

B. Moments of Z 



To be quantitative, let us choose a quantity M which is 
extensive meaning M = Nm where m is the "density" or 
per particle value. This is based on the additive property 
over subsystems M = ^ Mi . For a random system we better 
write M = M{N, {Q}), with {Q} representing all the random 
variables. To recover thermodynamics, we want [MJ^v to be 
proportional to A'^ for N ^ ex. Now, if it so happens that for 
large A*' 

M{N,{Q}) Nma, (16) 

with md independent of the explicit random variables, then 
M is said to have the self averaging property. Note that no 
averaging has been done in Eq. 1161 One way to guarantee 
this self-averaging is to have a probability distribution 

P{M/N) ^^-=5° S{md). (17) 



For thermodynamic quantities we know that it is the 
average of the free energy that is required. This is called 
quenched averaging. 

We may write 

[InZ],, = [ln{[Zl, + (Z - [ZU}1, 

= ln[ZU + I^%-^i^ + ... (20) 

to show the importance of the variance of the partition 
function. If the variance remains small, in the limit 
^ oo, then the polymer can be described by the av- 
erage partition function which is more or less like a pure 
problem. Otherwise not. 
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1. Hamiltonian for moments 

Higher order terms involve higher moments [^"]av oi' 
rather cumulants. These moments can be written as the 
partition function of an n-polymer problem with an extra 
interaction induced by the disorder. Starting from H as 
given by Eq. Q, and the gaussian distribution of Eq. 
©, we have 

n 

K]av = j ^''") ^(^)^^ 



[g-/3/i:'5(r.(^)-r),)(r,^)]^ 



exp (/32/AE,<, '^(r.(^) 




where iJ„ 



(3^Cdz E.<,'^(r.W-r,(z)) 



(21) 



This particular form can be understood in terms of two 
polymers. These two polymers start from the same point 
and do their random walk as they take further steps. 
If there is a site which is energetically favourable, both 
the polymers would like to be there. The effect is like 
an attractive interaction between the two polymers - an 
interaction induced by the randomness. 



2. Bound state: two polymer problem 

For the second moment, we have a two polymer prob- 
lem. The analogy with quantum mechanics tells us that 
for d < 2, any binding potential can form a bound state 
but a critical strength is required for a bound state for 
d > 2. In the polymer language, this means that any 
small disorder will change the behaviour of the free (pure) 
chain for d < 2 (disorder is always relevant) but for 
d > c?c = 2, if (3A < {(3A)c, the chain remains pure- 
like (disorder is irrelevant). Actually in higher dimen- 
sions (d > 2) the delta-function potential needs to be 
regularized appropriately (e.g., by a "spherical" well). 
Such cases are better treated by renormalization group 
(RG) which also helps in making the definitions of rele- 
vance/irrelevance more precise. We discuss this below. 

In short, the second term (fluctuation in partition func- 
tion) in the expansion of Eq. (|20|l cannot be ignored if 
d < 2 or if /3A is sufficiently large for d > 2. This sig- 
nals a disorder dominated phase for all disorders in low 



dimensions or at low temperatures (weak disorders) in 
higher dimensions. 



C. RG approach 

It is instructive to develop a renormalization group ap- 
proach for the two polymer problem. 



1. Expansion in potential 

We do an expansion in the interaction potential and 
just look at the first contributing term. Full series can of 
course be treated exactly. On averaging, the first order 
terms drop out, yielding 



dzP^Ad(ri(z) - r2(z)) + ...(22) 



This is the first order term if Eq. H21|) is used. A diagram- 
matic representation is often helpful in book-keeping as 
shown in Fig. Obut we avoid the details here. 



2. Reunion 

For this two polymer problem, the interaction con- 
tributes whenever there is a meeting or reunion of the 
two polymers at a site. At the order we are considering, 
there is only one reunion but this reunion can take place 
anywhere along the chain and anywhere in the transverse 
direction. 



FIG. 3: Renormalization of the interaction. The polymers 
are represented by the two lines and an intersection represents 
an interaction. 



The probability that two walkers starting from origin 
would meet at r at z is given by P^(r,z) (Eq. (0J) so 
that a reunion anywhere is given by a space integral of 
this probability which gives 



T^z= J drP^(r,z) 



(47rz)-*, with = d/2. (23) 



This exponent ^' will be called the reunion exponent. 
The occurrence of a power law is again to be noted. The 
eventual renormalization group approach hinges on this 
power law behaviour. 
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3. Divergences 

The contribution in Eq. H22() apart from some con- 
stants involve an integral over the reunion behaviour 
given in Eq. 1231 This integral in the limit N oo 
is 



N 



dz- 



1 



,d/2 



a}-'^!'^ for d > 2 (24) 
N^--^/^ for d<2. (25) 



For a finite cut-oflF, as is usually the case, the integral is 
finite for d > 2, and therefore [{Z - [Z]^^)\^ w 0{(3'^/^). 

This however is not the case for d <2 with d = 2 as a 
borderline case. The problem we face here is ideal for a 
renormalization group approach. 



4. RG flows 



d = 2). A marginally relevant variable introduces a new 
phase transition (or fixed point) in higher dimensions. 
The emergence of the new fixed point for d > 2 signals a 
phase transition. 



6. Fixed points: strong vs weak disorder 

Since u comes from the variance of the distribution, 
it cannot be negative. We therefore need to concentrate 
only on the u > part with the initial condition of u(L ~ 
a) = uq. What we see that for d < 2, the flow on the 
positive axis goes to infinity indicating a strong disorder 
phase for any amount of disorder provided we look at long 
enough length scale. An estimate of this length scale may 
be obtained from the nature of divergence for a given uq. 
An integration of the flow equation gives L ~ Uq^^^ '^^ 
(d < 2), a crossover length beyond which the effect of 
the disorder is appreciable. 



Let us introduce an arbitrary length scale L in the 
transverse direction and define a dimensionlcss "running" 
coupling constant 




u{L) ^ {f3^KA)L' 



(26) 



and study the RG flow of the coupling constant as the 
scale L is changed. This takes care of reunions at small 
scales to define the effective coupling on a longer scale 
(with rescaling). The flow equation is 



L— = {2-d)u + u^ 
dL 



(27) 



The magnitude of the coefficient of the term is not 
very crucial because at this order, this coefficient can be 
absorbed in the definition of the u itself. What matters 
is the sign of the term. 

Exercise III.l By change of variables in Eq. 1^21^ . show 
that u{L) as defined in Eq. \2b]) is the correct dimension- 
less variable. 



5. Relevant, irrelevant and marginal 

A coupling that grows (decays) with length scale is 
called a relevant (irrelevant) term because at long scales 
the contribution of this quantity cannot (can) be ignored. 
A coupling that does not change with scale is called a 
marginal variable. 

For example at the pure fixed point {u ~ 0), disorder 
is relevant (see Eq. (|26|) ) for d < 2 but becomes irrele- 
vant for d > 2, while it is marginal at d = 2. Marginal 
variables are important because the renormalization pro- 
cedure may add higher order corrections, like the term 
in Eq. H27|l . which then determines the growth or decay of 
the term. Here we see that the disorder is a marginally 
relevant variable, eventually increasing with length (at 



d = 1 d>2 

u= *^ u=o u^iir 



(a) 



(b) 



FIG. 4: RG Fixed points for u. (a) Based on the second 
moment of the partition function, (b) Based on the KPZ 
equation. Arrows show the flow of m. e = 2 — d. 

For d > 2, there is a fixed point at u* = |e| where 
e = 2 — d. For u < u* , the disorder strength goes to zero 
and one recovers a "pure"-like behaviour. This is a weak 
disorder limit. But, if u > u* the disorder is relevant. 



D. Ans.: Yes, it is! 

Based on the fixed point analysis, we conclude, as al- 
ready mentioned, that disorder can be relevant depend- 
ing on the dimensions we are in (i.e. the value of d) 
and temperature or strength of disorder. In particular, 
follows: 

1. A disorder-dominated or strong disorder phase for 
all temperatures for d < 2. 

2. A disorder dominated or strong disorder phase at 
low temperatures for d > 2. 

For d > 2, one sees a phase transition by changing the 
strength of the disorder or equivalently temperature for a 
given A. This is an example of a phase transition induced 
by disorder which cannot exist in a pure case. The fiow 
equation around the fixed point for d > 2 shows that one 
may define a "length-scale" associated with the critical 
point that diverges as ^ ~ |w— it*|~^ with C = c^l- In 

the weak disorder phase where the disorder is irrelevant, 
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[InZJjj^ w ln[Z']3^^. and therefore one may put a bound 
on the transition temperature Tc for a lattice model as 
Tc > Ta as defined below Eq. fl^ . 

In absence of any fixed point for the strong disorder 
phase in this approach, no further quantitatie results can 
be obtained about the phase itself. 



IV. QUANTITATIVE RESULTS ON 
EXPONENTS 



A strong disorder phase has characteristic sample fluc- 
tuations (described by the exponent 9). Even the size 
need not be gaussian like (i.e., ly ^ 1/2). We discuss 
various ways of getting these quantities, some of which 
are general and some are very particular for the directed 
polymer problem. 



A. Replica 

Even in the absence of any knowledge of the probability 
distribution like e.g. P{\nZ), various information can be 
extracted from the moments. This starts from the simple 
identity that 



in 



exp 



E 



1 



'[(InZ) 



Tni(c) 



(28) 

where [(lnZ)™]iv^ are the cumulants. As an interacting n 
polymer problem, each polymer of length iV, one expects 
an extensive term and a correction which depends on n 
and TV, as 



ln[Z"] 



nNe + T{n,N). 



(29) 



It is the correction term that is more useful and informa- 
tive. Our interest is generally in small n and — > oo, 
as shown in Fig. \^ This also becomes clear from a more 
familiar form for the free energy. Eq. [2H|can be written 
as 



lim 



1 



(30) 



so that to compute the average free energy we may con- 
sider a case of n-replicas of the original system or after 
averaging, an n-polymer problem with extra interactions 
induced by the disorder though an n limit is to be 
taken at the end. A few possible paths to take the limit 
for long chains are shown in Fig. [S| 



1/N 



computer /numerical (c) 
scaling path (b) 



Analytical approach (a) 



FIG. 5: Paths for replica approach 



Nontrivial results arc expected if and only if the origin 
in Fig. [SI is a singular point so that the limits n ^ and 
N ^ oo become non-interchangeable. In other words, 
the n and A^ dependences should be coupled so that the 
appropriate path is a scaling path like (c) in the figure 

Exercise IV. 1 Path dependence with rare events 

(i) Consider the probability distribution of Eqs. \2ci2b\ 
and calculate [Inx]^^. 

(ii) Calculate the moments [x^]^^ and use Eq. ^.WJ)) . Co 
along paths (a) and (c) of Fig. [3| and compare with the 
result of (i). 

(Hi) Show that correct average is obtained if a scaling 
path is chosen with u^/N = constant. 



1. 71 ^ and N oo 

Now, [(lnZ)^]iv is related to the free energy fluctua- 
tion and, as per Eq. l(Tl|l . it is expected to scale as N'^^ . 
If higher order fluctuations (or cumulants) do not require 
any new exponent, then it is fair to expect [(lnZ)™]iv'' ~ 
A^™^ A natural choice is T{n,N) = T{nN^). 

If we demand that In [Z"]^^ is proportional to A^ for 
large A^, then, for x = nN^ — > oo, J^{x) ~ x^^^ so that 



In [Z'\^ = neN + sn^/'^N. {N -> oo) 



(31) 



This is for path (a) of Fig. [5| In contrast, if we take 
n ^ for finite A^, path (c), a Taylor series expansion 
gives 



ln[Z"] 



nA^e + an^AT^^ + .... (n 0) (32) 



Eqs. I31I32I can be used to calculate 9, the free energy 
fluctuation exponent.. 

Exercise IV.2 Prove that In [Z"]^^ cx A^. 

The size of a polymer can also be handled in this replica 
approach by a careful limit as follows. 



[ < '^W > lav = 



/ VTZrlj exp(-/3H) 



= lim 



/x'7^cxp(-/3i^) 

(/ Vn exp{-/3H))"^^ J Vnr% exp(-/3i7) 



(33) 



Using the identity lim„^o [{/ Vnexp{-l3H)}'']^^ = 1, 
one may consider the partition function of the n-replica 
problem, compute the value of average of 7'^ (without 
normalization) for one of the replicas and then take n — > 
0. However this procedure is not easy to implement. 



B. Bethe ansatz 

The replica approach requires an evaluation of [Z'^]^^. 
For a gaussian distributed, delta-correlated disorder. 
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[Z^]g^ corresponds to the partition function of an n- 
polymer system with the Hamiltonian given by Eq. 
H21(l . Noting the similarity with the quantum Hamil- 
tonian with z playing the role of imaginary time, finding 
\n[Z"]^^ for — *■ oo is equivalent to finding the 
ground state energy i? of a quantum system of n par- 
ticles. This problem can be solved exactly only in one 
dimension (d = 1) using the Bethe ansatz. This gives 
the ground state energy as 



E = ~K{n-~r?) {d = \) 



(34) 



which gives 



{d=\) (see below) (35) 



The polymer has swollen far beyond the random walk 
or gaussian behaviour. What looks surprising in this ap- 
proach is that there is no "variance" (2nd cumulant) con- 
tribution. It is just not possible to have a probability dis- 
tribution whose variance vanishes identically! This is a 
conspiracy of the iV — > oo limit inherent in the quantum 
mapping and the value of the exponent B that suppressed 
the second cumulant contribution (see Eq. (|?T|l . To re- 
peat, in this quantum ground state approach, iV ^ oo 
has been taken first and so there are no higher order 
terms. 

Exercise IV. 3 Why cannot the Bethe ansatz he ex- 
tended to higher dimensions? 



1. Flory approach 

Using the quantum analogy, we may try to estimate the 
ground state energy in a simple minded calculation. The 
elastic energy is like the kinetic energy of quantum parti- 
cles which try to delocalize the polymers (random walk) 
while the attractive potential tries to keep the polymer 
together. For n polymers there are n{n — l)/2 interac- 
tions. We take the large n limit so that if the particles 
are bound in a region of size R, the energy is (using di- 
mensionally correct form with R as the only length scale) 



n 



(36) 



consistent with the 



n 

^"^m R 

which on minimization gives E ^ n 
Bethe ansatz solution. At this point we see the prob- 
lem of the replica approach if the limit is taken too 
soon. Since our interest is eventually in n ^ 0, we could 
have used in this argument the linear term of the com- 
binatorics. That would have made energy "extensive" 
with respect to the number of particle and replaced the 
disorder-induced attraction by a repulsion (note the neg- 
ative sign). The end result would however have no 
dependence. This is a real danger and any replica calcu- 
lation has to watch out of these pitfalls. Quite strangely 



we see that the correct answer came by taking n ^ oo 
first and then n — > or, probably better to say, by stay- 
ing along the "attractive part" of the interaction only. 

Exercise IV. 4 Does the Bethe ansatz require n ^ oo? 



V. CAN WE HANDLE FREE ENERGY? 

So far we have been looking at the moments of the 
partition function and trying to extract relevant infor- 
mation from that. Is it possible to avoid the replica trick 
altogether and treat the free energy explicitly? 

This is something unique to directed polymers that 
there is a way to study the average free energy and im- 
plement RG directly for the free energy bypassing the 
71 ^ problem of Fig|21completely. This is a very impor- 
tant step because it gives an independent way of checking 
the results of replica approach. 



A. Free energy and the KPZ equation 

For a polymer, the partition function satisfies a diffu- 
sion or Schrodinger-like equation. This equation can be 
transformed to an equation for the free energy i^(r, z) — 
— TlnZ(r, z). This is the free energy of a polymer whose 
end point at z is fixed at r. To maintain the distance fixed 
at r a force is required which is given by g = — Vi^. If we 
want to increase the length of a polymer by one unit, we 
need to release the constraint at the previous layer (think 
of a lattice). The change in free energy would then de- 
pend on the force at that point, and of course the random 
energy at the new occupied site. The change dF{r)/dz 
being a scalar can then depend only on the two scalars 
V • g, g^. A direct derivation of the differential equation 
for the free energy shows that these are the three terms 
required. The differential equation, now known as the 
Kardar-Parisi-Zhang equation, is 

^ = ^V^F-^(V^)^+.(r,.). (37) 

If we can solve this exact equation and average over the 
random energy rj^ we get all the results we want. 

One may also write down the equation for the force in 
this "fixed distance" ensemble as 

l! = §^'^-2]?^-^^+^^^^'^)- 

This equation is known as the Burgers equation. 



B. Thermal correlation and an exact result 

In the previous section we considered the force required 
to maintain the distance fixed. We now consider a fixed 
force ensemble that helps us in determining the response 
functions. 
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It is known in statistical mechanics that the response of 
a system in equihbrium is determined by the fluctuations. 
If we consider the response of a directed polymer to an 
applied force, the averaged response function comes from 
the Hamiltonian 



" - ¥ I - (I 



ri{v{z),z) -g - I —dz 



The disorder is gaussian-distributed as in Eq. 
The general response function for the force is 



dgidg. 



i,j representing the components. 



(39) 



[M-{n){r,)U, (40) 



C. Free energy of extension: pure like 

We want to know the free energy cost in pulling a poly- 
mer of length N from origin (where the other end is fixed) 
to a position r. For the pure case, the free energy follows 
from Eq. Q (with K inserted) as 

F(r, N) - F(o, N)^^K^. (pure) (44) 

For the disordered case, we use Eq. to see the 

surprising result, 

[Fir, N) - Fio, N)l, = ^ ^ ^- (disorder) (45) 

Therefore, on the average the stretching of a chain is 
pure-like (elastic) with the same elastic constant though 
the fluctuation is anomalous {6 7^ 0). This result has a 
far reaching consequence that in a renormalization group 
procedure, the elastic constant must remain an invari- 
ant. As we shall see, this invariance condition puts a 
constraint on 1/ and 6, making only one independent. 



Exact result on response: pure like 



D. RG of the KPZ equation 



By a redeflnition of the variables and using the S- 
correlation of the disorder in the z direction, wc have 



[lnZ(g)],, = [lnZ(g = 0)], 
from which it follows that 



IN 
d.K' 



TN 



(41) 



(42) 



as one would expect in a pure system, Eq. . And there 
are no higher order correlations. 

Two things played important roles in getting this pure- 
like result: (i) The disorder correlation has a statistical 
translational invariance coming from the delta function 
in the z-coordinate, and (ii) the quadratic nature of the 
hamiltonian. If disorder had any correlation along the 
length of the polymer, Eq. 021^111 not be valid. 

Exercise V.l Put a semi- flexibility term [d^r / dz'^Y in 
the Hamiltonian and then derive the expression for Ct ■ 

Exercise V.2 Is Eq. valid for a lattice problem? 

Note also that the full range [— cx),cx)] has been used. 



Exercise V.3 Prove that 

P{F{vm)) = P(F(0) 
Hint: Choose g = d.Kv/N . 



d.Ki 



N - 



2N 



(43) 



To analyze the nonlinear KPZ equation, an RG proce- 
dure may be adopted. This RG is based on treating the 
nonlinear term in an iterative manner by starting from 
the linear equation. This is a bit unusual because here 
we are not starting with a "gaussian" polymer problem, 
rather, a formal linear equation that does not satisfy the 
exponent relation Eq. (|48|) . Leaving aside such peculiar- 
ity, one may implement the coarse-graining of RG to see 
how the couplings change with length scale. 



1. Scale transformation and an important relation 

As a first step, we do the scale transformation 

x^bx,z^ b^l^z, and F b'^'^'F, (46) 
so that the randomness transforms like 

[77(ri,Zi)77(r2,Z2)]av 

= A(5(ri - r2)(5(zi - Z2) 

^ 6-'^-(i/-)A,5(ri - r2)(5(zi ~ Z2) (47) 

This transformation done on Eq. H45|l shows that for K 
to be an invariant (no 6-dependencc) we must have. 



1 = 1v 



(48) 



This is trivially valid for the gaussian pure polymer prob- 
lem but gives a relation between the free energy fluctu- 
ation and the size of the polymer. This is borne out 
by the intuitive picture we develop below. This relation 
gives the size exponent i^ = 2/3in(i = 1 (See Eq. (|35|) 'l. 
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The equation in terms of the transformed variables is 
then 

dz 2K 2K ^ ' 

^(i-d^-2e)/,.^(j.^ ^) (49) 

The 6-dependcnt factors can be absorbed to define new 
parameters, except for K. 

The temperature however gets rcnormahzed. Its flow 
is described by the flow equation 

L^— — T (to leading order) (50) 

For V > 1/2, T(L) 0. The disorder dominated phase 
is therefore equivalent to a zero temperature problem. 
In other words, the fluctuation in the ground state en- 
ergy and configurations dominate the behaviour at finite 
temperatures in situations with > 1/2. It is this renor- 
malization that was missing in Sec lIII Cl 

The nonlinearity (or g^) contributes further to the 
renormalization of the temperature through the appro- 
priate dimensionless variable u = [K A/T'^)L'^~''' (same 
as in Eq. H26() . The important flow equation is for this 
parameter u (upto constant factors) and it is 



2(i- 3 



L — = 2 - d)u . 



(51) 



One immediately sees a major difference with the flow 
equation Eq. H27|) for d = 1. Because of the change of 
sign of the quadratic term in Eq. 1511 there is now a f.p. 
for d = 1 as shown in Fig. 01b). This flow equation does 
not behave properly in a range d G [1.5,2) but that is 
more of a problem of RG than DP as such and so, may 
be ignored here. Note also that no extra information can 
be obtained for d > 2 other than what we have obtained 
so far in Sec. IIIII A point to stress is that the fixed point 
at d = 1 gives back the exact exponents we have derived 
so far. 



For the pure case as A'' ^ oo the width of P{r]y,N) 
increases. Hence the increase of Ct with TV. With ran- 
domness, for T —^ we need to look for the minimum en- 
ergy path. Let us suppose that there is a unique ground 
state, i.e. i5(rjv) or i^(rjv) is a minimum for a particular 
path. This tells us that as the temperature is changed T 
still low, the polymer explores the nearby region so that 
the probability distribution gains some width which is 
determined by the thermal length. Susceptibility would 
be the width of the distribution and this is independent 
of N. This cannot satisfy the relation. Since the relation 
is valid on the average, one way to satisfy it is to assume 
that most of the samples have a unique ground-state but 
once in a while (rare samples) there may be more than 
one ground state which happens to be far away from each 
other. Suppose there are such rare samples, whose prob- 
abilities decay as N^'^, where the paths are separated by 
N'^, then the contribution to the fluctuation from these 
samples would be N'^^~'^ . In case k = 2i/ — 1, we get back 
the exact result. The relation of Eq. H48|l tells vs n = 9. 
The rare events control the free energy fluctuation. 

What we see here that though the average behaviour 
is the same as that of the pure system, the underlying 
phenomenon is completely different; the average thermal 
response is determined by the rare samples that have 
widely separated degenerate ground state and the prob- 
ability of such states also decay as a power law of the 
chain length. 



In a given sample elastic energy 



The pinning energy would also grow with length say as 
. We sec, — 2v — 1 = 6. One way to say this is that 
the sample to sample fluctuation and the energy scale for 
a given sample are the same. 




FIG. 6: Various paths for various locations of the end point. 



VI. VIRTUAL REALITY: INTUITIVE PICTURE 

Powered by the quantitative estimates of the free en- 
ergy fluctuation and size exponents, we now try to gener- 
ate a physical picture. To do this we exploit a few more 
exact results. 



A. Rare events 



These results can now be combined for an image of 
the minimum energy paths. If the end points at z = N 
are separated by r ^ , the paths remain separated 
(each path exploring an independent disordered region) 
until they join at Az ~ r^^'^ ^ after which they follow the 
same path. If the end points are separated by a distance 
r ^ N'^ , the two paths explore independent regions and 
they need not meet. This (Fig|51) picture is often alluded 
to as the river-basin network. 



How to understand the exact result of Eq. (|^ ? We 
have shown that there is a low temperature region or 
in lower dimensions, disorder or randomness results in a 
new phase but then its response is the same as the pure 
system? 



B. Probability distribution 

For a pure polymer, the probability distribution of the 
end point is gaussian but it need not be so for the dis- 
ordered case. One way to explore the probability dis- 
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tribution is to study the response of the polymer as we 
take it out of its optimal or average position, e.g. by 
applying a force. In a previous section we used the fixed 
distance ensemble where the end point was kept fixed 
and we looked at the force g required to maintain that 
distance (see Eq. (|H8|I ) . Here we consider the conjugate 
fixed force ensemble. 



1. Response to a force 

Let us apply a force that tries to pull the end of the 
polymer beyond the equilibrium value r ~ i?o. In equi- 
librium, the average size R or extension by the force can 
be expressed in a scaling form 



R — Rn f 



.9^0 
T 



(52) 



This is because for zero force one should get back the 
unperturbed size while the force term may enter only in 
a dimensionless form in the above equation where the 
quantities available are the size Rq and the thermal en- 
ergy. For small g, linearity in g is expected. This require- 
ment gives 



R = Ro^g {kB = l), 



(53) 



R is not proportional to iV if ^ 1/2 (i?o ~ A^"). The 
polymer acts as a spring with T/Rq as the effective spring 
constant. 



2. Scaling approach 

Let us develop a physical picture and the correspond- 
ing algebraic description (called a scaling theory). The 
polymer in absence of any force has some shape of char- 
acteristic size Rq. The force stretches it in a way that it 
breaks up into blobs of size = T/g. For size < the 
polymer looks like a chain without any force and these 
blobs, connected linearly by geometry, act as a "new" 
polymer to respond to the force by aligning along it. We 
now have two scales Rq and ^g = T/g. A dimensionless 
form is then 



. S3 



T 



(54) 



Now each blob is of length Ng so that ^g ~ N'^ and there 
are N/Ng blobs. We therefore expect 



(55) 



This gives a susceptibility x = dR/dg ^ g 



(l-2i/)/i/ 



3. And so the probability distribution is 

Let us try to get the susceptibility of Eq. H55() in an- 
other way. Let us assume that the probability distribu- 
tion for large R is 

^5^ 



P(i?) ~exp(-(r/i?o) 



(56) 



The entropy is given by S{r) = — In P. The free energy 
in presence of a force which stretches the polymer to the 
tail region is given by 



F = Tir/RoY - gr. 



(57) 



This on minimization gives g ~ '^i^/Rf^Y ^ ■ By equat- 
ing this form with Eq. (|55|1 . we get 6 = 1/(1 — i/) and 



P{R) - exp 



r 

Rq 



(58) 



For 1/ = 1/2 we do get back the gaussian distribution. 




FIG. 7: A blob picture of the polymer under a force. Though 
drawn as sphere, the 2-direction is elongated with isotropy in 
the transverse direction. 

The above analysis, done routinely for polymers, re- 
lies on the fact that there is only one length scale in the 
problem, namely, the size of the polymer. If we are en- 
titled to do the same for the disorder problem, namely 
only one scale, i?o ~ N'^, matters, then the blob picture 
goes through in toto. The chain breaks up into "blobs" 
and the blobs align as dictated by the force. Each blob is 
independent and the polymer inside a blob is exploring 
its environment like a DP pinned at one end. The proba- 
bility distribution is therefore given by Eq. 1581 which for 
d = 1 is 

P{r) ^ cM~\r\Vt') {d=l). (59) 

If we use the relation AF ee F{x, N) - F{0, N) - x^/N, 
then the above probability distribution can be mapped 
to the distribution of the free energy as 

1/2(1-!^)" 



P(F) 



exp 



exp 



\AF\ 

NO 

AP|3/2 

7VV3 



(d = l). 



(60) 
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C. Confinement energy 

Suppose we confine the polymer in a tube of diameter 
D. This is hke the locahzation length/ argument used 
to justify the energy in the quantum formulation. The 
polymer in the random medium won't feel the wall until 
its size is comparable to that, D Nq which gives the 
length at which the polymer feels the wall. Elastic energy 
of a blob is D'^/Nq. But because of the tube, the polymer 
will be stretched in the tube direction. The number of 
blobs is N/No, so that the energy is 

This gives the known form used in the quantum 

analogy (consistent with dimensional analysis) but for 
1/ = 2/3, this gives 1/D. 

A cross-check of this comes from the energy of a blob. 
Each blob has the fluctuation energy Nq and so free en- 
ergy per unit length F/N - N^/Nq - i:)-2(i-i')/i^. 

VII. OVERLAP 

The problem we face in a disordered system is that 
there is no well defined ground state - the ground state 
is sample dependent. There is therefore no predefined 
external field that will force the system into its ground 
state (e.g., a magnetic field in a ferromagnetic problem). 
This is a generic problem for any random system. 

But, suppose, we put in an extra fictitious (ghost) poly- 
mer and let it choose the best path. Now we put in the 
actual polymer in the same random medium but with a 
weak attraction v with the ghost. At T = 0, this DP will 
then sit on top of the ghost. In absence of any interac- 
tion {v ^ ), the DP would go over the ghost in any 
case if there is a unique ground state, not otherwise. At 
non-zero temperatures there will be high energy excur- 
sions and how close to the ground state we are will be 
determined by the number of common points of the two 
polymers. This is called overlap which may be quantita- 
tively defined as 

g. = ^y'<5(ri(z)-r2(z)), (62) 

for a given sample i and then one has to average over the 
disorder samples, q = [gijav 

In case of a repulsive interaction, the situation will be 
different. If there are more than one ground state, the 
two chains will occupy two different paths and there is no 
energetic incentive to collapse on top of each other when 
the repulsion v 0+. In such a scenario, the overlap 
q{v 0+) 7^ q{v 0—). Conversely, a situation like 
this for the overlap would indicate presence of degenerate 
ground states. For a unique ground state, the second 
chain would follow a nearby excited path with certain 
amount of overlaps with the ground state. 



In a replica approach, such occupancy of different 
ground states may be achieved by "replica symmetry 
breaking" (i.e. various replicas occupying various states) 
but the difficulty arises from the n —^ limit. In the 
case of directed polymer, we have argued that the degen- 
erate states occur only rarely and therefore the effect of 
"replica symmetry breaking" if any has to be very small. 
This is why the Bcthe ansatz gave correct results without 
any breaking of replica symmetry. 

The method to compute the overlap was developed 
by Mukherji. By introducing a repulsive potential 
V J dzS{ri{z) — r2(z)) for the two polymers in the same 
random medium, the free energy F(ri(z), r2(2:), z, v) can 
be shown to satisfy a modified KPZ type equation 

f - E - J^i^.Ff + ,(r. .)) +.^(rr-r.). 

i ^ ^ 

(63) 

which can be studied by RG. The mutual interaction has 
no effect on the single chain behaviour but the interac- 
tion gets renormalized. The fiow equation for the dimen- 
sionless parameter u of Eq. 1511 remains the same. The 
exponent relation of Eq. 1481 also remains valid. The in- 
teraction gets renormalized as 

dv (1-6 , u\ 

where v is in a dimcnsionless form. For the pure prob- 
lem {6 = {],v = 1/2) this reduces to the expected flow 
equation of Eq. (|27|l for repulsive interaction {u —> —u). 
For overlap one needs only the first order term because 
we need f — > 0. 

The overlap can be written in a polymer-type scaling 
form q^ f{vN-'t"'), where E = 6'-0i/-l. The above 
RG equation for v shows that the exponent E = at the 
stable fixed point for u of Fig^Jb) at d = 1. However, 
E < at the transition point for d > 2. This means 
that the overlap vanishes at the transition point from 
the strong disorder side as g ~ |T — Td'^l'' . 

This approach to overlap can be extended to m-chain 
overlaps also, which show a nonlinear dependence on m 
at the transition point. This suggests that eventhough 
the size exponent is = 1/2 gaussian like, there is more 
intricate structure than the pure gaussian chain. 

VIII. SUMMARY AND OPEN PROBLEMS 

The behaviour of a directed polymer in a random 
medium in 1-1-1 dimensions seems to be well under- 
stood. There is a strong disordered phase at all temper- 
atures. For d > 2 renormalization group analysis shows 
a phase transition from a low temperature strong disor- 
dered phase to a weak disorder, pure-like phase. There 
are rare configurations with degenerate widely separated 
ground states, giving a contribution to "overlap". How- 
ever most questions on the higher dimensional strong- 
disorder phase remain open. 
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